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1.  INTRODUCTION 

Providing  environmental  support  for  Navy  operations 
requires  special  attention  to  the  marine  boundary  layer.  Ac¬ 
curate  meteorological  nowcasts  and  forecasts  play  a  key  role 
in  weapons  system  support,  aircraft  operations,  and  optimum- 
track  ship  routing.  Furthermore,  surface  wind  stresses  and 
heat  fluxes  from  the  atmospheric  model  are  used  to  drive 
ocean  thermodynamic  and  hydrodynamic  models  (Clancy. 
19S7).  Through  these  coupled  air/ocean  systems,  the  atmo¬ 
sphere  indirectly  plays  a  role  in  search  and  rescue  missions 
and  submarine  operations. 

The  Navy  Operational  Global  Atmospheric  Prediction 
System  (NCGAPS)  consists  of  four  separate  components — 
data  quality  control,  data  analysis,  initialization,  and  the 
forecast  model.  NOGAPS  operates  in  a  6-hr  update  data 
assimilation  cycle,  with  the  6-hr  model  forecast  serving  as 
the  background  field  for  the  data  analysis,  which  in  turn  pro¬ 
vides  the  fieids  for  the  initialization  of  the  foreevt  model. 
The  forecast  model  is  an  18  7 -level.  79-wave  triangular  trun¬ 
cation.  global  spectral  model  with  a  nonlinear  normal-mode 
initialization  (Hcgan  and  Rosmond.  19911.  The  data  analysis 
uses  the  multivaiate  optimum  interpolation  (MVOI)  tech¬ 
nique  described  cy  Goerss  and  Phoebus  (1991)  to  analyze 
geopotential  height  and  the  u-  and  e-wind  component  incre¬ 
ments  on  standard  oressure  surfaces.  The  analysis  fields  arc 
produced  for  16  standard  pressure  levels  {1006-10  mb.  plus 
925  mb)  on  the  1.5"  resolution  Gaussian  grid  of  the  spectral 
forecast  model. 

Since  the  MVOI  is  a  three-dimensional  incremental  anal¬ 
ysis.  ob  •valions  made  at  various  heights  can  influence  the 
analysis  results  at  1000  mb.  However,  high  quality  upper- 
air  profile  measurements  over  the  oceans  are  limited  to  a  few 
island  radiosonde  stations  and  an  occasional  balloon  launch 
from  a  ship.  While  these  observations  provide  some  infor¬ 
mation  about  the  structure  of  the  lower  atmosphere,  their 
numbers  are  too  few  to  have  much  impact.  The  conven¬ 
tional  upper-air  observations  are  supplemented  by  temper¬ 
ature  soundings  from  the  polar-orbiting  satellites.  Although 
these  soundings  are  plentiful,  accuracy  is  a  problem.  Other 
data  sources  include  single-level  wind  reports  from  aircraft 
and  cloud-tradced  wind  observations  estimated  from  the  geo¬ 
stationary  satellites.  The  low-level  satellite  winds  are  nor¬ 
mally  concentrated  at  the  925-830  mb  level,  while  the  upper- 
level  satellite  winds  and  the  aircraft  data  arc  predommenily 
at  the  250-200  mb  level.  Over  water,  conventional  sources 
of  surface  data  include  observations  from  ships,  moored  and 
drifting  buoys,  and  Coastal- Marine  Automated  Network  sta¬ 
tions.  These  marine  surface  reports  are  concentrated  primar¬ 
ily  in  the  northern  hemisphere,  along  major  shipping  routes. 

The  relative  sparseness  of  data  over  the  oceans  certainly 
hinders  efforts  to  produce  accurate  analyses  and  forecasts  of 
meteorological  conditions  in  these  areas.  The  effect  is  most 
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noticeable  in  the  southern  hemisphere,  of  course,  where  the- 
ratio  of  sea  to  land  is  greater  and  the  lack  of  data  is  even  more 
pronounced.  The  only  practical  solution  to  this  problem  is 
to  develop  and  deploy  high-quality  remote  sensors.  A  step 
towards  that  goal  was  achieved  on  June  19.  1987,  when  the 
first  Special  Sensor  Microwave /Imager  (SSM/I'i  was  launcned 
by  the  Defense  Meteorological  Satellite  Program  iDMSP). 
The  SSM/I  is  a  passive  microwave  radiometer  that  provides 
high-resoiution  measurements  of  a  variety  of  environmental 
parameters — among  them,  ram  rate,  precipitable  water,  ma¬ 
rine  wind  speed,  sea  ice  location  and  concentration,  sou  mois¬ 
ture,  and  land  surface  temperature. 

Of  particular  interest  for  meteorological  da’s  assimila¬ 
tion  are  the  SSM/I  wind  speed  data.  While  the  resolution, 
quality,  and  coverage  of  the  SSM/1  data  are  excellent,  their 
usefulness  is  somew  hat  limited  by  the  lack  of  directional  infor¬ 
mation.  Nonetheless,  early  tests  indicated  that  there  would 
be  a  small  positive  impact  if  the  data  were  included  it,  the 
analysis  (Goerss.  19S9a).  Given  that  this  positive  impact 
would  be  in  an  area  of  high  Navy  interest,  nameiv  the  ma¬ 
rine  boundary  layer,  the  operational  assimilation  of  SSM/1 
wmd  speed  data  into  NOGAPS  was  initiated  on  Septemoer 
12.  1990.  The  positive  impact  of  these  observations  on  the  op¬ 
erational  system  is  described  in  Phoebus  and  Goerss  ;1991). 
The  initial  implementation  of  the  SSM/I  wind  speed  date 
used  a  ,-reiiminary  1000  mb  analysis  to  assign  wind  direction. 
Since  that  lime,  we  have  investigated  alternative  methods  of 
assigning  wind  direction,  and  those  results  are  reported  here. 

2.  SSM/I  DATA  PROCESSING 


The  SSM/I  data  are  processed  and  stored  by  the  Navy- 
in  near  real-time  at  Fleet  Numerical  Oceanography  Center 
(FNOC).  FNOC  uses  the  algorithm  of  Goodberiet.  et  of. 
(1990)  to  estimate  the  ocean  surface  (19.5  ml  wind  speeds 
from  the  SSM/I  measured  brightness  temperatures.  These 
speeds  have  a  range  of  3-25  m  s'1 ,  and  meet  the  specified 
accuracy  of  ±2  m  s"1  in  rain-free  conditions.  Since  heavy- 
concentrations  of  water  vapor  and  rain  effectively  mask  mi¬ 
crowave  emission  from  the  ocean  surface  at  the  SSM/I  fre¬ 
quencies,  wind  speed  retrievals  under  these  conditions  are 
more  likely  in  error.  Currently,  the  NOGAPS  data  pre¬ 
processor  screens  the  SSM/I  data  so  that  only  rain-free  obser¬ 
vations  are  used  in  the  MVOI.  Furthermore,  any  wind  speed 
reported  outside  the  range  of  the  SSM/I  instrument  is  ig¬ 
nored. 

The  wind  speed  data  are  also  subjected  to  quality  con¬ 
trol  procedures  that  compare  the  observed  wind  speeds  to  a 
preliminary  analysis  field.  Whether  or  not  an  observation  is 
rejected  during  this  comparison  depends  upon  noth  the  ob¬ 
served  wind  speed  and  the  wind  speed  from  the  preiiminar- 
analysis  valid  at  the  observation  location.  The  strictest  cr.te- 
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ria  are  applied  and  the  most  data  are  rejected  5%)  when 
the  analyzed  wind  speed  is  very  light  and  the  SSM/1  wind 
speed  is  substantially  stronger.  Since  the  preliminary  analy¬ 
sis  will  also  be  used  to  assign  wind  directions  to  the  SSM/I 
observations,  this  decision  reflects  the  fact  that  it  would  be 
inappropriate  to  ai><gn  ;«  -.ji.d  Hir»**  tjon  from  a  light  arid 
variable  condition  to  a  stronger  wind  speed.  The  oilier  cate¬ 
gories  exhibit  more  tolerance  for  differences  between  the  pre¬ 
liminary  analysis  and  the  SSM/1  data  and.  as  a  result,  reject 
fewer  observations. 

The  preliminary  analysis  is  run  using  the  same  software 
and  the  same  background  fields  a<  the  NOGAPS  analysis, 
but  is  executed  at  only  2  hours  pas:  the  synoptic  time  ?  com¬ 
pared  to  -*-9-5  hours  for  the  update  cycle),  so  it  has  a  shorter 
data  window.  Analyzed  fields  are  produced  only  at  the  four 
pressure  levels  nearest  the  surface.  The  preliminary  analysis 
uses  almost  all  of  the  available  data,  except  for  aircraft  wind 
reports,  satellite  temperature  soundings,  and.  of  course,  the 
SSM/I  wind  speed  data. 

2.2  Computation  of  Sunerobs 

A  typical  G-hour  coverage  by  the  SSM/I  produces  in  ex¬ 
cess  of  200.000  observations  at  approximately  25  km  resolu¬ 
tion.  a  resolution  much  higher  than, the  effective  resolution  of 
NOG  APS-  Thus,  to  reduce  the  number  of  data  points  without 
losing  information  on  scales  relevant  to  NOGAPS,  the  SSM/I 
data  are  combined  ~ i.U.  u.e  •mtiaiiy  local 

averages  of  groups  of  observations. 

Before  the  superobs  are  computed,  the  wind  speeds 
from  the  preliminary  analysis  are  interpolated  to  the  loca¬ 
tion  of  each  individual  observation,  and  the  interpolated  val¬ 
ues  are  subtracted  from  the  observed  wind  speeds  to  form 
SSM/I  wind  speed  increments.  Next,  the  globe  is  divided 
mto  200  km  boxes  and.  given  sufficient  numbers  of  obser¬ 
vations.  the  SSM/I  wind  speed  increments  within  each  box 
are  averaged.  The  increments  rather  than  tne  original  val¬ 
ues  are  used  because,  in  general,  the  increment  field  is  mor* 
horizontally  homogeneous  than  the  raw  data.  Once  tne  su- 
perob  increments  are  computed,  the  superob  wind  speeds  are 
obtained  by  adding  the  increment  in  each  box  to  the  prelim¬ 
inary  analysis  wind  speed  interpolated  to  the  center  of  the 
box. 

3.  ASSIGNMENT  OF  WIND  DIRECTION 

The  analysis  variables  are  geopotential  height  and  the 
u-  and  t>wind  components — variables  that  have  well-known 
dynamical  relationships.  The  design  of  the  MVOI  analy¬ 
sis  is  such  that  it  cannot  utilize  scalar  wind  speed  informa¬ 
tion  directly.  So.  once  the  SSM/I  superob  wind  speeds  have 
been  computed,  each  superob  wind  speed  must  be  assigned 
a  wind  direction.  Tb^re  have  been  several  techniques  pro¬ 
posed  for  independently  determining  appropriate  directions 
for  the  SSM/I  wind  speeds.  Atlas  ct  al.  1 1991)  combine 
conventional  wind  data  and  SSM/I  winds  with  background 
analyzed  surface  wind  fields  from  the  European  Centre  for 
Medium  Range  Weather  Forecasts,  using  variational  tech¬ 
niques  to  constrajn  the  degree  to  which  the  SSM/I  data  can 
modify  the  background.  Yu  ( 1987)  bases  his  technique  on  Ek- 
man  boundary  layer  dynamics,  using  the  National  Meteoro¬ 
logical  Center's  surface  pressure  analysis  and  the  SSM/I  wind 
speed  to  compute  the  surface  drag  coefficient,  which  then  al¬ 
lows  the  computation  of  the  vector  wi~d  at  some  height  z* 
wuhjn  the  boundary  layer  {typically  10  m). 

In  our  case,  we  use  the  power  of  a  global  data  assim¬ 
ilation  system  to  produce  a  preliminary  analysis  using  the 


most  recent  data,  independent  of  the  SSM/I.  Originally,  the 
1000  mb  wind  field  from  the  preliminary  analysis  was  used  to 
assign  SSM/I  wind  direction.  However,  the  NOGAPS  fore¬ 
cast  model  produces  a  10  m  wind  field  from  the  iowest  model 
ct -level  (-N.  40  m).  primarily  to  provide  forcing  data  for  the 
FNOC  oceanographic  models.  The  10  m  winds  are  sub;ht> 
dependent  and  consistent  with  the  model's  oounaarv  layer 
physics.  Thus,  the  assigned  SSM/I  wind  directions  are  now 
taken  from  the  analyzed  10  m  wind  field,  which  is  horizontally 
interpolated  to  the  location  of  each  superob  Once  the  direc¬ 
tions  have  been  assigned,  the  SSM/I  superobs  are  processed 
like  any  other  surface  wind  observation  and  are  assumed  to 
nave  the  same  observational  error  as  wind  observations  from 
ships. 

An  improved  technique  for  assigning  wind  direction  to 
the  SSM/I  wmd  speeds  is  under  development.  The  cur¬ 
rent  operational  technique  just  described  is  heavily  dependent 
upon  the  MVOI  analysis  background  fields  m  areas  sparsely 
covered  by  marine  surface  data.  The  accuracy  of  these  back¬ 
ground  fields,  which  are  G-hour  forecasts  from  the  global  spec¬ 
tral  model,  is  dependent  upon  the  performance  of  the  forecast 
model.  In  the  northern  mid-latitudes,  where  weather  systems 
typically  move  from  west  to  east,  the  accuracy  of  the  oack- 
ground  fields  over  ocean  areas  is  greatest  just  to  the  east  of 
the  data-rich  continental  areas  and  decreases  from  west  to 
east.  The  largest  errors  m  the  background  fields  are  associ¬ 
ated  *!iii  incorrect  fuie».A>vi  oi  tne  position  and/or  intensity 
of  extra- tropical  cyclones.  The  error  that  is  most  critical  to 
the  assignment  of  SSM/I  wind  direction  is  that  associated 
with  an  incorrect  position  forecast  In  such  situations,  tne 
background  fields  normally  depict  quite  well  the  structure  of 
the  cyclone  and  its  associated  fronts  but  snow  an  offset  :n 
their  positions.  As  a  result,  in  the  absence  of  other  data, 
the  wind  directions  assigned  to  the  SSM/I  wmd  speed  oQser- 
vatiotis  by  the  current  technique  would  reflect  the  erroneous 
offset.  Ti  c  goal  of  our  research  is  to  utilize  the  SSM/I  wmd 
speed  observations  to  determine  and  correct  for  such  offsets 
before  the  wind  direction  is  assigned. 

Prior  to  computing  cross- soectrai  quantities  like  the  co¬ 
herency  or  phase  spectra,  it  is  common  practice  in  time  se¬ 
nes  analysis  to  align  the  two  time  senes  (Jenkins  and  Watts. 
1968).  This  is  done  by  lagging  one  time  senes  with  respect 
to  the  other  and  computing  the  values  of  the  cross- covariance 
function  for  the  different  lags.  The  lag  at  which  the  maxi¬ 
mum  value  of  the  cross-covariance  function  is  found  is  used  to 
align  the  time  series.  This  technique  was  utilized  by  Duchon 
and  Goerss  ( 1976)  in  their  time  series  analysis  of  aircraft  data 
from  intercomparison  Sights  performed  during  the  National 
Hail  Research  Experiment,  These  flights  typically  consisted 
of  three  aircraft  flying  in  formation  with  one  in  the  lead. 
The  data  from  these  flights  consisted  of  samples  1  second 
apart.  Before  performing  cross-spectrum  analysis  upon  the 
time  series  from  the  aircraft,  the  data  set:  fr?n  *.b*  trailing 
aircraft  were  aligned  to  that  from  the  lead  aircraft  by  mea¬ 
suring  the  lime  displacement  from  zero  of  the  peak  m  the 
cross- covariance  functions.  The  location  of  this  peak  could 
be  easily  defined  and  the  resulting  time  displacement  was 
attributed  to  the  persistent  longitudinal  separation  between 
the  aircraft  and  to  ^synchronization  of  the  aircraft  docks. 

We  propose  to  use  an  analogous  technique  in  areas 
sparsely  covered  by  marine  surface  observations  to  spatially 
align  the  SSM/I  wind  speed  observations  with  the  prelim¬ 
inary  analysis  wind  fields  prior  to  the  assignment  of  wmd 
direction.  For  an  area  of  appropriate  size,  for  which  there 
»s  sufficient  coverage  by  SSM/I  wmd  sp^ed  observations,  we 
can  compute  the  correlation  between  the  wmd  speed  observa- 
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uons  and  the  wind  speeds  from  the  preliminary  10  m  analysis 
at  the  location  of  the  observations.  We  then  spatially  shift 
the  preliminary  analysis  with  respect  to  the  SSM/1  observa¬ 
tions  and  recompute  the  correlations  for  spatial  lags  in  the 
east-west  and  north-south  directions  to  determine  the  lag  for 
which  this  correlation  is  a  maximum.  The  wind  directions 
from  the  shifted  analysis  are  then  assigned  to  the  SSM/1  wind 
speed  observations,  but  the  observations  themselves  are  not 
actually  displaced. 

We  illustrate  this  technique  for  an  area  located  between 
30° S  and  60°S.  30®  W  and  120°W.  In  Fig  1  the  isotachs  from 
the  unshifted  preliminary  10  m  analysis  are  plotted  along 
with  the  SSM/1  wind  speed  suberobs.  In  the  northern  half  of 
this  domain,  there  is  a  northwest-to-southeast  axis  of  mini¬ 
mum  wind  speeds  in  the  background  field  that  is  located  to 
the  southwest  of  a  similar  axis  in  the  wind  speed  observations. 
Similarly,  the  center  of  maximum  winds  in  the  south-centra: 
part  of  the  domain  is  located  southwest  of  the  corresponding 
maximum  wind  speed  superobs.  Table  I  provides  the  corre¬ 
lation?  computed  by  systematically  shifting  me  preliminary 
analysis  relative  to  the  SSM/1  wind  observations,  using  1.5® 
increments  in  the  east-west  direction  and  1.0®  increments  in 
the  north-south  direction.  The  maximum  correlation  is  found 
by  shifting  the  preliminary  analysis  1.0®  to  the  north  and  4.5" 
to  the  east. 


Fig  I.  SSM/1  wind  speed  superobs  fits)  available  at  12Z.  August 
22,  1991 ,  compared  to  the  preliminary  analysis  of  10  m  wind  speeds 
coun  toured  every  j  kts. 


TABLE  1.  Cross-correlations 
Between  SiM/I  and  Offset  Background  Wind  Speeds 
Com du ted  over  Area  from  60*-30®S.  120* -90* W 
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Fig.  2  shows  the  analyzed  1000  mb  wind  field  computed 
using  the  current  technique,  with  direction  assigned  from  tc* 
preliminary  10  m  wind  analysis.  In  comparison.  Fig  3  suo*. 
the  comparable  field,  but  in  this  case,  the  SSM/1  wind 
uons  were  assigned  by  shifting  the  background  held  reia->\* 
to  the  data.  As  you  can  see  from  the  streamline*.  tfc:«  r*  • 
m  a  slight  eastward  shift  m  the  centers  of  circulation,  mo,- 
noticeable  in  the  northwesterly  section  of  the  end 


Fig  2  Streamlines  and  isotachs  .contoured  every  10  kts j  at  luOC 
mo.  from  analyzed  fields  produced  for  J2Z.  August  22.  1991  um* 
SSM/1  wind  data  mih  directions  assigned  from  the  preliminary  1C 
m  analysts  field 


Fig.  3  Same  as  Fig  2.  but  with  urn d  direction  assigned  to  the 
SSM/1  data  from  the  shifted  10  m  analysis  field  most  highly  corre¬ 
lated  with  the  SSM/1  data  over  (he  area  mown. 


We  plan  to  test  the  application  of  this  technique  nrsi 
m  the  southern  hemisphere  mid -latitudes.  Th*  marine  sur¬ 
face  data  coverage  in  the  northern  hemispnere  mid-rAtituce? 
is  fairly  extensive  and  is  supplemented  bv  rhe  generation  o. 
synthetic  surface  observations  m  the  v»cimtv  ol  exira-tropica. 
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cyclones  (Goerss.  1989b).  While  the  tropics  are  not  as  well 
covered  by  marine  surface  data,  significant  tropical  circula¬ 
tions  arc  accurately  depicted  in  NOGAPS  through  the  inclu¬ 
sion  of  synthetic  soundings  that  are  generated  using  informa¬ 
tion  contained  in  the  warning  messages  issued  by  the  Joint 
Typhoon  Warning  Center,  Guam  and  the  National  Hurri¬ 
cane  Center,  Miami  (Goerss  et  al.,  1991).  Thus,  the  largest 
errors  are  undoubtedly  associated  with  extra-tropical  highs 
and  lows  in  the  Southern  Hemisphere.  However,  the  tech¬ 
nique  may  prove  useful  in  other  areas,  as  well. 

While  this  experiment  demonstrates  the  readability  of 
the  technique,  for  operational  application  we  must  determine 
ways  to  automatically  identify  those  areas  where  a  phase  er¬ 
ror  in  the  background  is  likely.  In  that  regard,  this  work 
is  closety  related  to  another  of  our  research  interests,  that  is. 
identifying  areas  where  the  background  errors  are  larger  than 
those  assumed  by  the  MV01  analysis,  so  that  the  specified 
background  error  can  be  modified  in  real-time  over  limited 
areas.  The  goal  of  this  effort  is  to  identify  quantities  that 
serve  as  useful  predictors  for  locating  such  occurrences.  Once 
these  areas  are  Identified,  they  would  also  be  likely  candidates 
w  here  the  alternative  method  of  assigning  SSM/I  wind  direc¬ 
tion  would  prove  beneficial. 

■1.  CONCLUSIONS 

SSM/T  wind  speed  data  provide  an  abundant  source 
of  information  for  improving  low-levet  marine  analyses  and 
forecasts,  an  area  of  key  importance  for  support  of  naval  ot>- 
erations  and  the  ocean  modeling  community.  All  indications 
are  that  this  data  source  is  of  high  quality  and  has  had  a 
positive  impact  on  the  performance  of  the  global  data  as¬ 
similation  system.  To  utilize  the  SSM/I  data,  the  primary 
decision  'hat  must  be  made  is  how  to  assign  a  wind  direction 
to  each  observation.  In  NOGAPS,  we  currently  use  the  10  m 
wi.-.d  fields  computed  from  ?  preliminary  low-ievei  analysis  of 
conventional  data  and  satellite  doud-lrackcd  winds. 

However,  in  spite  of  the  relatively  close  agreement  be- 
tween  the  SSM/T  data  and  the  preliminary  analysis  wind 
speeds,  visual  comparison  of  the  fields  indicate  that  in  some 
cases  there  are  phase  differences.  A  preliminary  study  com¬ 
puting  cross-correiaiions  of  the  SSM/1  and  preliminary  anal¬ 
ysis  wind  speeds  substantiates  this  phase  error.  When  the 
SSM/I  data  'locations  are  offset  in  another  direction,  the 
cross-correlation  with  the  background  is  often  increased,  while 
the  number  of  SSM/I  observations  rejected  by  the  analysis 
quality  control  procedures  is  reduced.  Thus,  future  efforts 
will  concentrate  on  improved  ways  to  assign  wind  direction 
to  the  SSM/I  wind  speeds  and  how  to  implement  such  proce¬ 
dures  operationally.  These  methods  should  also  prove  useful 
for  resolving  any  wind  direction  ambiguity  problems  associ¬ 
ated  with  the  ERS-l  scatterometer  winds  when  they  become 
available. 
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